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Abstract. High- multiplicity proton-proton collisions at the LHC may exhibit 
collective phenomena such as elliptic flow. We study this issue using DIPSY, a brand- 
new Monte Carlo event generator which features almost-NLO BFKL dynamics and 
' describes the transverse shape of the proton including all fluctuations. We predict the 

. eccentricity of the collision as a function of the multiplicity and estimate the magnitude 

of elliptic flow. We suggest that flow can be signaled by a sign change in the four- 
particle azimuthal correlation. 



The observation of high-multiplicity events in pp collisions at the LHC opens up 
an interesting possibility that the collective flow, usually discussed in the context of 
nucleus collisions, may be realized in the final state of pp collisions [lj. Indeed, highest 



multiplicity events in the 7 TeV pp run have dN c h/dr] ~ 40 which is comparable to 
semi-central Cu-Cu collisions at RHIC, and flow has been observed in the latter. In 
this contribution we study some key questions about the possibility to observe elliptic 
flow j2] in pp using DIPSY — a recently released Monte Carlo event generator [3]. More 
details can be found in [lj. For related works, see, [3 El [7]. 

Firstly, one should bear in mind that high-multiplicity pp collisions and nucleus 
collisions with the same multiplicity are vastly different. In particular, the former are 
rare fluctuations in the broad N C h distribution while the latter refer to average events 
at fixed centrality. There are several sources of multiplicity fluctuations in pp: (i) 
Impact parameter fluctuation — Unlike in nucleus collisions, in pp collisions it is not 
possible to determine the impact parameter for each event. High multiplicity events 
mostly come from collisions with small impact parameter, (ii) Intrinsic fluctuation of 
the proton's wavefunction — Protons at the LHC undergo the QCD evolution to become 
a dense system of small-x gluons. Since the QCD evolution is stochastic, there are 
large event-by-event fluctuations in the gluon number. High multiplicity events arise 



Eccentricity and elliptic flow in pp collisions at the LHC 



2 



from protons with an unusually large occupation number of gluons. (iii) Fluctuation 
in the collision process — In a single high-multiplicity pp event, there are many (more 
than 10) gluon-gluon scatterings. The number of subcollisions fluctuates due to the 
probabilistic nature of collisions (partonic cross section), (iv) Fluctuation in the final 
state parton showering — High multiplicity events typically contain several jets, and the 
fragmentation of jets is a stochastic process. 

At first, point (i) seems to be a fatal blow to any hope of observing elliptic 
flow t>2 in pp. Naively, in central collisions the eccentricity would be very small, 
hence small, unobservable v^- However, this may be solved by point (ii). The QCD 
evolution generates fluctuations not only in the gluon number, but also in the transverse 
distribution of gluons because the gluon splitting probability depends on the transverse 
coordinates!]] This makes it possible to have a sizable participant eccentricity even in 
central collisions 



where a XjV are the variances of x and y coordinates of liberated gluons. [By 'liberated 
gluons' we mean gluons which actually interacted as well as those in the underlying 
events.] We have computed p]) using DIPSY [3J [8], a new event generator which takes 
into account all the points (i)-(iv) above. It is based on the QCD dipole model [9] which 
is the coordinate space formulation of the BFKL evolution, and therefore captures the 
correct transverse dynamics of small-x gluons. In addition to the leading-order BFKL, 
DIPSY features a dominant part of the next-to-leading corrections, energy conservation, 
and saturation effects. Thanks to this, the energy dependence of observables is a 
prediction in DIPSY. Once the parameters have been fixed at one value of energy, 
there is no ad hoc re-tuning of parameters at different energies. 

The result for the eccentricity e{2} = yj (ejLrt) and e{4} = (2( e 2 art ) 2 - (e 4 ,^)) 1 / 4 
as well as the interaction area (S) at 7 TeV is plotted in Fig. [T] (left) as a function of 
N c h within the ALICE acceptance \rj\ < 0.9. [(...) denotes event-by-event averaging in 
a given multiplicity bin.] We see that the eccentricity is about 30-40% in the highest 
multiplicity region. This is similar to the value in semi-central nucleus collisions at 
RHIC and the LHC. 

Next we give an estimate t> 2 . In nucleus collisions, v% and e are roughly proportional 



where X/Kq is a certain parameter fitted to experimental data [10]. This empirical 
formula works both at RHIC and at the LHC, and both for Au-Au and Cu-Cu at 
RHIC although they differ in size by a factor of two. The latter supports the general 
argument that the applicability of hydrodynamics is controlled by the dimensionless 

X Previous works considered the fluctuation of 'hot spots' [6] and 'flux tubes' [7]. There the transverse 
distribution of these objects was assumed to be random. In our case the transverse distribution of 
gluons is not random, but governed by the QCD evolution. 
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Figure 1. Predictions for the eccentricity and the interaction area S (x0.1)[fm 2 ] 
(left) and elliptic flow (right) versus the charged multiplicity in the interval |?7| < 0.9 
at y/s = 7 TeV. 



parameter a = dN/drj ra ^ er than the system size, and the necessary condition a < 1 
is well satisfied in high-multiplicity pp events. We thus expect that (121) gives a reasonable 
estimate of t>2 even in pp collisions. The result for t>2{2} and t>2{4} are plotted in Fig. [1] 
(right). [In this calculation we actually used a slightly improved version of (j2j), see, jl].] 
Elliptic flow is about 6%, comparable to the value found at the LHC. 

Is it possible to observe the flow contribution v 2 ~ 6%? Experimentally, v 2(2, 4} 
are measured from the azimuthal angle correlation 

v 2 2 {2} = (cos(2(0 4 - 0,))) , (3) 
v 4 2 {4} = 2(v 2 {2}) 4 - (cos(2(0, + ^ - <f> k - m . (4) 

They differ from the genuine v-i by the so-called nonflow contribution (v2{n}) n = v% + 5 n 
where 5 n is the n-particle correlations not associated with flow. In nucleus-nucleus 
collisions, they are relatively innocuous because they scale with the multiplicity as 
S n ~ l/N^ 1 . In pp collisions, however, one expects significant nonflow contributions 
from various initial and final state effects, and the ALICE collaboration has indeed found 
a much slower decrease in N c h [H]. This slow decrease is also observed in Monte-Carlo 
simulations, meaning that the isolation of the flow contribution v% could be difficult 
for small values of n. In the highest multiplicity bins, ^2(2} ~ 0.13 |TT] which is 
twice as large as the flow contribution. This implies that the two-particle correlation is 
dominated by nonflow effects. 

We are thus led to turn to higher order cumulants V2{n} with n > 4 which are 
by definition insensitive to two-particle nonflow correlations. Actually, the ALICE 
collaboration has also measured t>2{4} [H]. It turns out that the magnitude of t>2{4} 
in the data is still larger than our flow prediction, but very interestingly, it has the 
'wrong' sign — the rhs of (jlj) is negative! The same phenomena can be seen in PYTHIA 
and DIPSY (without flow effects), see Fig. [2j On the other hand, in the flow scenario 
(v 2 {4}) 4 is positive. [Note, however, that the statistical error bars on the preliminary 
measured (t>2{4}) 4 [U] are an order of magnitude larger than the value we predict 
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Figure 2. The DIPSY result for the fourth cumulant of the nonflow correlation 
QC{4} = -(w 2 {4}) 4 compared with a rough parametrization of the preliminary ALICE 
data p]. 

from flow.] We thus propose to look for this sign change in experiment as a possible 
signature of flow: If there is flow in the large N c h region, then the fourth order cumulant 
01]), which is negative in the absence of flow, will decrease in magnitude and eventually 
turn positive. 
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